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O Abstract 

l^-j We study the electroweak phase transition by lattice simulations of an effective 

<3\ 3-dimensional theory, for a Higgs mass of about 35GeV. In the broken symmetry 

phase our results on masses and the Higgs condensate are consistent with 2-loop 
perturbative results. However, we find a non-perturbative lowering of the transition 
O i temperature, similar to the one previously found at m# = SOGeV. For the sym- 

metric phase, bound state masses and the static force are determined and compared 
with results for pure SU (2) theory. 



There are strong indications that the electroweak standard theory predicts a first order 
phase transition at the electroweak scale [l]-[8]. The generation of the baryon asymmetry 
of the universe at the electroweak phase transition is an exciting hypothesis. A better 
quantitative characterization of the electroweak phase transition is required, however, in 
order to clarify whether the known baryon asymmetry could actually be generated in this 
way. 
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The electroweak phase transition cannot be treated completely by perturbative tech- 
niques. Within this frame, the symmetric phase would consist of massless W bosons. 
Then there are infrared problems which prevent a perturbative evaluation of the free en- 
ergy to higher loop order. The true, non-perturbative behaviour of the symmetric phase 
will be characterized by massive W- and Higgs bound states instead of massless W gauge 
bosons, however. Condensate formation may also lower the free energy of the symmetric 
phase. Therefore, non-perturbative techniques are required in order to determine the crit- 
ical temperature at which the free energies of both phases are equal. The free energy of 
the broken symmetry phase can be evaluated perturbatively, as long as the Higgs mass is 
not too large. If Higgs particles are heavy the broken phase will also be strongly coupled, 
for temperatures close to the critical one. This problem will not be considered in the 
present paper. Our analysis is limited to a study of the phase transition for m# in the 
range of 3bGeV. Our aim is to shed light on the non-perturbative behaviour of the sym- 
metric phase, and to test the (dimensionally reduced) method through a comparison with 
the perturbative treatment of the broken phase where the latter should be appropriate. 

The electroweak phase transition has been investigated by 3 + 1-dimensional lattice 
simulations before [4, 5, 7]. Since the interesting dynamics is expected to be carried by the 
static (constant in r) low momentum modes, dimensional reduction should be possible, 
also in the neighbourhood of the phase transition. In this case lattice simulations could 
be done with an effective 3-dimensional action generated by integrating out the non-static 
modes including the fermionic ones. Other than in the case of QCD, dimensional reduction 
should work for the electroweak theory in the interesting temperature range because g 2 is 
small. For that reason, the non-local 3-dimensional effective action may be approximated 
locally, even near to the critical temperature. In connection with the electroweak phase 
transition this approach has been pioneered by Farakos et al. [3, 6, 8]. 

We study the 5C/(2)-Higgs system with one complex doublet of variable length. The 
gauge field is represented by the unitary 2x2 link matrices U x ^ and the Higgs fields are 
written as & x = p x V x (p x = \Tr{Q+Q x ) is the scalar 'Higgs length', V x an element of the 
group SU(2)). The lattice action is 

S = fa 5> - \TrU p ) -PhY, \Tr{^ x U x ,^ x+ ,) + Y(pI + M ~ !) 2 ) (!) 



with 



Pa = — 2 - (2) 



In three dimensions the lattice Higgs self-coupling (5 R is given by 

A3 P% 
9l Pg 



Pr = ^1T- ( 3 ) 
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g\ and A 3 denote the 3 — d continuum gauge and Higgs self couplings which are 3 — d 
renormalization group invariants. They are related to the corresponding four dimensional 
couplings via 

g! = g 2 T, \ 3 = T(\ + 0(g*)) . (4) 



The string operators of length L = 1, 2, ... are defined as 

El(L) = ^t^x^Ux+ix^Ux+2ix,ix---U x+ (L-l)ix,ix^x+Lfj, (5) 

and used to form Higgs and VF-operators. Actually the 3 — d masses (inverse correlation 
lengths) have been obtained from the connected correlators between separated equal 'time' 
(in 2 + 1 dimensions) planes of the sum of 'spatial' string operators within these planes 
projecting out the proper SU(2) and spin content. To increase the overlap with the lowest 
mass states we choose L = 4 for the VF-mass determination and also partly in the Higgs 
mass sector (there also pi is used). As expected, the signal is cleaner and stronger in the 
low temperature phase. To determine the static force we have used Wilson loops W(R, T) 
of asymmetric extensions 5 < R < T < N/2 (N 3 is the lattice size). 

The vectorized Monte Carlo algorithm combines a three-dimensional Gaussian heat 
bath for the gauge fields and a four-dimensional Gaussian heat bath for the Higgs field. 
To reduce the autocorrelations near to the phase transition a heat bath step was followed 
by several reflections (eight in practice) for the Higgs and one reflection for the gauge 
field. 

We consider the quartic coupling = 0.0239. According to the tree level based 



relation 



A 3 A 1 m 2 H 



~2 ri '2 Q nm 2 



(6) 



#3 9 ° m w 

this would correspond to the case of zero temperature masses m# = 3bGeV, mw = 
80GeV. One should notice, however, that eq. (6) is only approximately true. It depends 
on details of the dimensional reduction (loop order, corrections for the adjoint Higgs field 
which we neglect altogether). Comparison with other calculations should therefore be 
done at the corresponding A value, and not necessarily refering to m# = 3bGeV. The 
lattice Higgs self coupling (5 R is not fixed in this setting but runs with (3 H according to 
eqs. (3,6). 

In the simulations we used 20 3 and partly 30 3 lattices. Typically 20,000 to 100,000 
iterations at the smaller lattices and up to 30,000 iterations at the larger ones have been 
performed for every coupling pair ((3h,(3r) at given j5o (with correlator measurements 
every 10 iterations). Near to the phase transition the largest integrated autocorrelations 
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have been observed in observables containing Higgs field degrees of freedom. We were 
cautious to have the measured correlation lengths (up to 10) on the 20 3 lattices not to be 
strongly influenced by the lattice size. We checked a few of the largest correlation lengths 
on larger lattices and found no significant lattice size behaviour. A particular example is 
presented later in Fig. 2 . 

Note that we never have to specify the value of the 4 — d gauge coupling g 2 . Masses 
are obtained in units of g 2 T after eliminating the lattice constant a through eqs. (2,4). 

We now turn to a discussion of results. Near to a first-order phase transition the two 
phases will coexist. If the transition is strong enough the two phases should manifest 
themselves as a two-state signal in the histograms of suitable observables. In Fig. 1, we 
show a histogram for < p 2 > (which is the length squared of the Higgs field p 2 averaged 
over the lattice for each measured configuration). For the particular lattice parameters 
chosen (/5 G = 16, flu = 0.33930) the peak of the broken phase has obviously more weight 
indicating that we are still slightly below the critical temperature. Locating the critical 
temperature precisely is a difficult problem. As it is known from previous studies [4] 
different indicators like the specific heat or the Binder cumulant result in various pseudo- 
critical couplings in a finite volume. A detailed finite size analysis is necessary to locate 
the inifinite volume critical parameters and to avoid these ambiguities. 
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From the position of the peak of the broken phase we can determine the renormalized 
value of the Higgs condensate, according to 
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«p 2 > - (pIJ) 



(7) 



where {p 2 ym ) is the expectation value of p 2 in the symmetric phase (close to the critical 
point) where it is actually almost independent of /5#. The relation (7) neglects a small 
temperature dependent two-loop contribution. In subtracting p 2 of the symmetric peak 
we define the condensate to be zero in the symmetric phase. This cannot be strictly 
maintained because vr is not renormalization group invariant and depends on the 3 — d 
renormalization scale p^ at two-loop order. The resulting condensate value is very small 
for any reasonable choice of p 3 , however. 

Equation (7) allows us to measure directly the quantity ^ which appears as coupling 
in quasiclassical estimates of the sphaleron rate. The sphaleron rate is sufficiently small 
(in order to avoid washing out any generated asymmetry) as long as ^ is of the order of 
| (or smaller) at the critical temperature. 

The best way to compare the predictions of continuum perturbation theory for the 
phase transition in any detail with lattice results would be to calculate the effective 
potential to two loop order, using both ways of regularization. This is not yet available. 
Therefore, we compare the relation between two observables in the two approaches. This 
does not require to relate explicitely the bare mass squared of the continuum formulation 
tttJ to the gauge-Higgs coupling (3 H . We also do not have to introduce the temperature 
explicitely. 

Fig. 2 gives a summary of our mass measurements at (5g = 12. With the operators 
under consideration we find two massive bound states in the symmetric phase, instead of 
the massless W bosons one would naively expect from perturbation theory. As long as the 
system tunnels between the two phases a safe measurement of the Higgs correlation length 
is not possible while the lightest IT^-state is not influenced by this tunneling. By inspecting 
the Monte Carlo histories of various operators it was possible to extract subsamples to 
perform correlation measurements which refer to the pure phases. This was the way to 
clearly identify a jump of the Higgs mass between the low and high temperature phases, 
respectively. In the IT^-channnel we were not able to establish an analogous discontinuity 
within our errors. Measurements on the larger lattice (30 3 ) do not indicate severe finite 
size effects near to the transition. There is also no indication that a heavy state becomes 
much lighter going to a larger volume. This would be the way a light state would appear 
in a small volume expansion. 
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Higgs and W masses in lattice units vs. gauge Higgs coupling (5h 



Table 1: Measured VF-masses at critical /?# vs. Pc on a 20 3 lattice 



Pg 


Ph 


m w a 


m w /(g 2 T c ) 


12 


0.34138(2) 


0.482(16) 


1.446(48) 


16 


0.33930(4) 


0.333(17) 


1.332(68) 


20 


0.33806(4) 


0.261(14) 


1.305(70) 



To identify finite a effects one has to compare masses at the same temperature. To get 
an idea about the strength of these effects we compare measured VF-states at couplings 
Ph which are, within our precision, nearest to the critical ones and, therefore correspond 
to the same critical temperature. For each /5 G , we select the suitable Ph value according 
to the largest correlation length in the Higgs sector. The data shown in Table 1 are 
consistent with each other for (5c = 16 and 20 while that at (5c = 12 significant finite a 
corrections cannot be ruled out. One should bear in mind, however, that in the present 
stage we cannot make sure that all these measurements correspond to exactly the same 
temperature with sufficient accuracy. 

We now discuss the mass measurements in the broken phase. One important quantity 
to judge the reliability of perturbation theory is the wave function renormalization for the 
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Lattice results for the wave function renormalization of the W field vs. (gT)/vR compared 
to one-loop perturbative predictions 



W field. It can be measured by relating the W mass to the renormalized condensate vr 

m w = Z w 2 ^Y^. (8) 

Results are shown in Fig. 3, together with the one-loop prediction calculated in 
Feynman gauge [9]. Within the statistical errors agreement is found. Z w is close enough 
to one to expect that the broken phase itself can be understood perturbatively in this 
range of light Higgs masses. Again, the data points refering to (3q = 12 are systematically 
somewhat below the perturbative curve, which can be interpreted as indication for a finite 
a correction at the strongest gauge coupling that we have investigated. 

The Higgs mass constitutes the other piece of information we can use to test the 
perturbative behaviour of the broken phase. Results are shown in Fig. 4, for — plotted 
versus the corresponding measured value of Within errors, we again find agreement 
with perturbative results obtained from the Feynman gauge two-loop potential [10]. The 
perturbative Higgs mass is obtained through the second derivative of the potential at the 
broken minimum (i.e. at zero momentum). This might account for a slight but systematic 
trend of the perturbative mass ratio to be larger than the lattice data which correspond 
to pole masses. The perturbative masses are renormalized by the corresponding one-loop 
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Measured mass ratio as function of (gT)/vR compared to perturbation theory 



wave function renormalizations Z h 2 and Z w 2 , resp., in Feynman gauge taken from refs. 
[9,10]. 

Finally, we discuss the static force. Results are shown in Fig. 5 for two values of 
0h corresponding to the broken and symmetric phase, respectively. Using the measured 
Wilson loops to determine the static potential and the force from Creutz ratios would 
result in an overestimation of the errors. Therefore we have fitted the large 'time' depen- 
dence of prolongated Wilson loops by an exponential in order to determine the ground 
state energy of the pair. Using these fitted energies with their errors we have defined the 
static force. The measured static force obtained from the naive Creutz ratios are in good 
agreement with this prescription. In Fig. 5, the perturbative one-loop contribution in the 
broken phase 

is also shown. The measured value of m w has been plugged into this expression. For the 
3 — d gauge coupling, two values have been used: the bare g\ as well as an effective coupling 
g\ R describing best the data. This effective coupling sums up finite loop contributions 
and is in good agreement with the reported wave function renormalization Z w . 
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Fig. 5 

Static force vs. distance in symmetric (open circles) and broken (filled circles) phase at 
0G = 12. The solid (dotted) line is the perturbative one loop contribution with g\a = 
0.41(0.33), m w a = 0.485. 



Whereas the static force is well represented perturbatively in the broken phase, a 
non-perturbative contribution is clearly seen in the symmetric phase. Roughly, it can be 
described as a linear contribution to the potential, with a string tension being in good 
agreement with a = 0.11(g 2 T) 2 found by Teper (0.13(# 2 T) 2 at f3 G = 12) [11] in 3 - d 
pure SU(2) gauge theory. This confining-type behaviour should not extend to very large 
distances beyond a screening length. Screening is not seen, however, up to the distances 
we could measure. 

The remaining feature to be discussed in connection with the phase transition is an 
eventual lowering of the critical temperature due to nonperturbative effects (condensates) 
present in the symmetric phase. Returning to Fig. 4 we notice the perturbative two-loop 
estimate of the phase transition to occur at about ^ = 0.46. This is the endpoint of the 

line. Lattice results point towards ^ = 0.41 instead. We can now simply calculate the 
free energy that corresponds to that latter value to 2-loop accuracy in the broken phase. 
This value can then be identified with the non-perturbative (condensate) contribution to 
the free energy of the symmetric phase. The result is 

AF = -0.027(# 2 T) 3 (10) 
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somewhat depending on the gauge (results shown represent an average of Landau and 
Feynman gauge). This value is very sensitive to the value of ^ considered to be the 

critical one. ^ = 0.42 would give AF = — 0.019((/ 2 T) 3 instead. This should be compared 
with the value of AF = — 0.016((/ 2 T) 3 given in ref. [6] for m H = SOGeV. Thus, we confirm 
the presence of such a term which seems to be more or less independent of the Higgs self- 
coupling. Therefore, it should essentially be understood in terms of the dynamics of 3 — d 
pure gauge theory alone. 

In summary, our lattice results are in reasonable agreement with perturbative pre- 
dictions for the Higgs phase, as they should. The symmetric phase shows a confining 
behaviour at the intermediate distances we have explored. The string tension is very 
close to the value found for the 3 — d pure SU(2) theory. Two massive bound states are 
identified in the symmetric phase. The mass of Higgs state is smaller than that of the 
lightest glueball of the pure SU (2) theory. The lowering of the critical temperature, pre- 
sumably induced by gauge condensates in the symmetric phase, shows little dependence 
on the Higgs self-coupling A. 
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